Late stage kinetics for various wicking and spreading problems 
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The kinetics of spreading of a liquid drop in a wedge or V-shaped groove, in a network of such 
grooves, and on a hydrophilic strip, is re-examined. The length of a droplet of volume f2 spreading 
in a wedge after a time t is predicted to scale as fi 1 / 5 ^ 5 , and the height profile is predicted to be 
a parabola in the distance along the wedge. If the droplet is spreading radially in a sparse network 
of V-shaped grooves on a surface, the radius is predicted to scale as f2 1//6 i 1//3 , provided the liquid is 
completely contained within the grooves. A number of other results are also obtained. 

PACS numbers: 68.08.Bc, 47.10. +g, 05.45.-a 
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I. INTRODUCTION 

Wetting in complex geometries and on rough surfaces 
provides a wealth of fascinating non-linear hydrodynam- 
ics problems, as well as being of commercial importance 
in numerous industrial sectors. Perhaps the first kind of 
problem to be considered was the penetration of liquid 
into porous materials, where Washburn in 1921 demon- 
strated that the distance attained by the wetting front 
follows a i 1 / 2 law where t is time [l|. Much later the 
spreading of droplets on flat surfaces was addressed by 
various workers, such as Tanner and Lopez et al al- 
though it took some time for the subtleties of the physics 
at the wetting front to be resolved 0, 0, 0, ■ Gener- 
ally, the wetting front advances with a t a law where a 
is a small exponent which depends on the geometry of 
spreading and the origin of the driving force. For exam- 
ple, a = 1/10 for a drop spreading radially driven by 
surface tension (Tanner's law), and a — 1/8 for a droplet 
spreading radially driven by gravity (see Oron et al Q 
for a summary of results). 

The kinetics of wetting on rough surfaces has also been 
investigated experimentally and theoretically H, 0, 0] . 
A paradigm for this problem is the spreading of a liquid 
in a wedge or V-shaped groove 0, 0, ^| ; indeed wet- 
ting in a network of V-shaped grooves has been invoked 
recently for oil spreading on skin Another kind of 

problem that has been considered is the wetting of hy- 
drophilic strips 0| i as an example of wetting in a con- 
trolled microstructure that might be contemplated in a 
microfluidic device. In all these problems, a t 1 ! 2 spread- 
ing law has been observed, but in the cases considered 
thus far, there has been a reservoir which provides liquid 
at essentially a constant pressure. In the present paper, 
the problems of spreading in a wedge, in a network of V- 
shaped grooves, and on a hydrophilic strip are revisited. 
It is found that in the absence of a reservoir, the spread- 
ing law changes to t a with a < 1/2, similar to Tanner's 
law and related problems. 

These problems are first approached by scaling argu- 
ments developed in the next section. The bona fides of 
the scaling arguments is established by rederiving some 
known results for spreading on flat surfaces. In a further 
section, the scaling exponents are recovered by similar- 



ity analysis on the underlying partial differential equa- 
tions which govern spreading. This also allows the scaling 
shape of the spreading drops to be computed. 



II. SCALING ARGUMENTS 

The Washburn problem of a liquid being drawn into a 
capillary tube of internal dimension d shown in Fig.^a) 
is considered first pj . This is a model for penetration of 
liquid into a porous material for which d interpreted as 
a mean pore size. The arguments here are very famil- 
iar, but form the basis for the more complex problems 
considered below. 

Once the liquid has penetrated a sufficient distance 
L ^S> d, a Poiseuille law obtains for the liquid velocity 
and the penetration rate, thus 



d 2 Ap 

T] L 

where 77 is viscosity, and the pressure drop 



dL 
~dt 



Ap~ a/d 



(1) 



(2) 



is due to the surface tension a of the curved surface at 
the wetting front, at a mean curvature ~ 1/d. All ge- 
ometric factors associated with the shape of the tube 
and a finite contact angle have been dropped, although a 
contact angle 8 < tt/2 is required for imbibition to take 
place. Combining Eqs. JIJ and J5J gives 



dL 
~dt 



a d 



(3) 



which integrates to 



(atd/rj) 1 / 2 . 



(4) 

This is the simplest form of the Washburn equation pj . 
The result arises from a constant pressure drop acting 
over an increasing length of liquid, which responds by 
flowing according to the Poiseuille law. As we shall see 
below, this can be the case for many situations where a 
reservoir of liquid is present, but if a reservoir is absent, 
the rate of spreading can be much slower. 
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Next, the problem of a drop of liquid spreading on a flat 
surface is considered, as shown in Fig. ^b) . Usually the 
problem is approached by an appeal to the hydrodynam- 
ics in the vicinity of the moving contact line 0, , but 
it can be analysed using similar concepts to the Wash- 
burn problem. Whilst only previously known results are 
recovered, the approach serves to illustrate further the 
arguments that will be used for the other problems. 

Consider a drop of liquid spreading on a flat surface, in 
the case of complete wetting. Let a measure of the radius 
of the spreading drop be R and the height in the centre be 
h. In the lubrication approximation, assuming a scaling 
shape of the droplet, all velocities will be proportional to 
(h 2 /rj)(Ap/ R) (compare Poiseuille law above) where Ap 
is the pressure drop between the centre and the radius 
R. In particular the drop radius is expanding at a rate 



dR 
~dt 



h 2 Ap 



(5) 



First consider the capillary spreading case where the 
pressure gradient is due to surface tension a. Simple 
geometry shows that the mean curvature at the centre of 
the droplet for h <C R is ~ h/ R 2 therefore the pressure 
drop is 



Ap ~ ah/R 2 



and hence 



dR 
~dt 



r\ R 3 ' 



(6) 



(7) 



If h were to be constant, as in the Washburn problem, 
this would be enough to determine the spreading rate. 
Here, though, a second relation connecting R and h must 
be sought. The desired relation follows from the conser- 
vation of drop volume f2, 



hR 2 ~ n. 



Thus 



which integrates to 



dR 
~dt 



R-iattf/rj) 1 ' 10 . 



(8) 



(9) 



(10) 



This result is Tanner's law The basic scaling R ~ 
573/10^1/10 j g we vj documented and has been experimen- 
tally verified 0. 

For the case where the spreading is driven by gravity, 
one has Ap ~ pgh where p is the mass density and g 
is the acceleration due to gravity. Following the same 
line of argument as above, one obtains R ~ (pgtfl 3 /!-]) 1 ^ 8 
. The behaviour crosses over from capillary spreading 
to gravity spreading when the Bond number pgR 2 /a in- 
creases. Since R is increasing, this means that capillary 
spreading always crosses over to gravity spreading if one 



$ d 



(a) 



L 



00 



2R 





w 



FIG. 1: Various wetting problems: (a) wicking into a cap- 
illary, (b) spreading on a flat substrate, (c) spreading in a 
wedge, and (d) spreading along a hydrophilic strip. 



waits long enough. The weak increase in spreading rate 
has been observed experimentally |8(- 

Another case that can be considered is planar or one- 
dimensional spreading. The only thing which changes 
is the volume conservation law which becomes hL ~ f2 
where L replaces R as the measure of extent of spreading, 
and is a volume per unit length. This yields L ~ 
(crffi 3 /??) 1 / 7 and L ~ {pgtSfi/rj) 1 / 5 for capillary @] and 
gravity |3( spreading respectively. 

In the next problem, exactly analogous arguments are 
applied to the case of spreading in a wedge, shown in 
Fig. Die). I n t ne case of spreading from a reservoir, 
this problem has been addressed by Romero and Yost 
[l2j |. The basic idea is that one has scale invariance, 
with the depth h of fluid being the only relevant length 
scale. Hence the transverse curvature of the interface 
oc 1/h. Thus, provided the droplet has become suffi- 
ciently extended so that the contribution of the longuitu- 
dinal curvature to the mean curvature can be neglected, 
the pressure p oc {—)a/h where the negative sign obtains 
if the surface is convex into the liquid. This is the case if 
26+<f> < 7r where 9 is the contact angle and 4> is the wedge 
angle as in Fig. ^c). In this case, the pressure becomes 
more negative as the amount of fluid is in the wedge gets 
smaller. This provides a pressure gradient which drives 
the liquid from regions of high loading to low loading. 

Even though the liquid has a free surface, a Poiseuille- 
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like law obtains, 



dL 
~dt 



h 2 Ap 
7~ 



(11) 



where L is a measure of the extent of spread of the liquid 
drop. The pressure drop follows from the above argu- 
ments, Ap ~ cr/h, and therefore 



dL 
~dt 



a h 



(12) 



in perfect analogy to Eq. J2J in the Washburn problem. 
The difference here is that h is a dynamic variable which, 
similar to the Tanner's law derivation above, is found 
by a volume conservation law. In this case AL ~ 
where A ~ h 2 is the cross section area occupied by the 
liquid, and is the volume of liquid. Thus h ~ (fl/L) 1 ' 2 . 
Substituting into Eq. i|12|) results in 



which integrates to 



dL a il 1 ' 2 
~dt ~ rj Ul 2 



(am 1 ' 2 /!!) 2 ' 5 . 



(13) 



(14) 



This spreading law, L ~ fl 1 ' 5 ^' 5 , is a new prediction. 
It can also be recovered by a more detailed analysis of 
the underlying partial diffential equations which will be 
given in a later section. 

Note that if the liquid had been spreading along the 
wedge from a reservoir, this would correspond to h ~ 
constant. This leads to the Washburn result L ~ t 1 ' 2 
and is the origin of the scaling law obtained by previous 
workers [HG3. 

The next case to be considered is the problem of a 
droplet spreading into a network of grooves. This has also 
been considered by various groups [a, 0, . If spreading 
occurs from a reservoir, then the front advances with a 
t 1 ' 2 Washburn-like law. However the case where the liq- 
uid is completely confined in the grooves is different. The 
scaling law in this case follows from arguments similar to 
those already applied to a drop spreading in a wedge. 
The analysis assumes that the grooves are rather sparse 
on the surface, in particular that the volume of liquid 
in the junction zones can be neglected compared to the 
volume contained in the grooves. 

Consider therefore a drop of liquid spreading in a 
sparse random network of grooves. It will spread essen- 
tially radially. Suppose that the groove cross section is 
a V-shape, so that both the capillary pressure and the 
Poiseuille scaling laws can be taken over from the case 
of spreading in a wedge. The fact that the grooves are 
randomly inclined with respect to the radial pressure gra- 
dient only introduces an additional numerical prefactor 
Thus Eq. (O above still holds (with L replaced 
by R). What changes is the volume conservation law: as 
R increases, more grooves become filled. If the length 



of grooves per unit area is Z -1 , where I is a character- 
istic groove spacing on the surface, the total length of 
grooves occupied by the liquid ~ R 2 /l and the volume 
fl ~ h 2 R 2 /l. Eliminating h between this and Eq. ifT^I) 
(with L replaced by R) gives dR/dt ~ (a/^ilfl) 1 ' 2 / R 2 . 
This integrates to 



R-iatim) 1 ' 2 /^' 3 . 



(15) 



The basic prediction therefore is that the spreading rate 
should slow from the t 1 ' 2 Washburn-like law for spread- 
ing from a reservoir, to an R ~ Q 1 ' 6 ^' 3 law as the liquid 
becomes confined to the grooves. The result is confirmed 
by a more detailed analysis of the underlying partial dif- 
ferential equations given later. 

The final problem that is considered is spreading along 
a hydrophilic strip, shown in Fig. ^d). The case where 
spreading occurs from a reservoir has been investigated 
both theoretically and experimentally, and it is found 
that the spreading front advances with a Washburn-like 
t 1 ' 2 law [l5j|. The situation in the absence of a reservoir 
was not investigated though. Once again, the absence of 
a reservoir leads to a slower rate of spreading, and the 
scaling law can be determined using analogous arguments 
to those applied above. 

Consider a drop of liquid which is completely wetting 
on a hydrophilic strip of width w. In the late stages, 
let the height of the liquid above the surface be h <C 
w, and a measure of the extent of spreading be L > 
w. The capillary pressure which drives the spreading is 
due to the transverse curvature of the interface, which 
for h <C w is ~ h/w 2 . The pressure gradient is thus 
Ap/ L ~ ah/w 2 L. In the lubrication approximation with 
h <C w, a Poiseuille- like law obtains for the fluid velocity 
with h being the relevant length scale. Hence the rate of 
extension of the droplet obeys dL/dt ~ (h 2 /rj)(Ap/ L) ~ 
(a / r))h 3 / w 2 L . Volume conservation indicates hwL ~ il 
and eliminating h between this and the above spreading 
rate finds dL/dt ~ (a / 'ry)Sl 3 '/ 'w 5 L 4 . This integrates to 
L ~ (atQ 3 /^ 5 ) 1 ' 5 . This is the predicted scaling law 
for the late stages of spreading in this problem, in the 
absence of a reservoir. 



III. SIMILARITY METHODS 

The results obtained above can also be derived by us- 
ing similarity methods to analyse the underlying partial 
differential equations (see for example Ref. |f|). Focus 
first on the problem of spreading in a wedge shown in 
Fig. nic). An equation which expresses local conserva- 
tion of liquid in the wedge is 



dA 
~dt 



d{Av) 
dx 







(16) 



where A ~ h 2 is the local cross section area occupied by 
liquid, v the mean velocity of the liquid, and x is distance 
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FIG. 2: Shapes of droplets spreading in (a) a wedge and (b) a 
network of grooves. These are similarity solutions g(u) found 
by integrating Eq. \Y2'2t with boundary conditions g = A and 
g' = at u = 0, for A = 1(1)5, where u is a scaled distance. 
In case (a) the shape is a parabola given by Eq. I|25|l . In case 
(b) the shape is obtained by numerical integration. 



along the wedge. The Poiseuille law indicates that the 
mean velocity follows 



v oc 



h 2 dp 
rj dx 



The arguments above show that p oc —a/h thus 

dp a dh 
dx h 2 dx 



(17) 



(18) 



Combining Eqs. I|16|) - (|18|) gives the following equation 
for the time evolution of the depth of liquid in the groove 
(compare Eqs. (8a)-(8c) of Romero and Yost [12j) 



d(h 2 



= K- 



d 



,dh\ 



•q dx V dx) 



(19) 



The dimensionless coefficient K (6, (f>) is given by Eq. (8c) 
in Romero and Yost 

El 

in terms of the static contact 
angle 8 and the wedge angle <j). For spreading to occur, 
K > is required. This corresponds to 29 + <f> < ir, or 
a liquid interface which is convex into the liquid. For 
the remainder of the discussion, the factor Ka/r/ and 
other trivial numerical prefactors can be adsorbed into 
the units of time, and will be omitted. 



Eq. (|19(l is basically a non-linear diffusion equation and 
one can seek similarity solutions of the form 



h(x,t) ~ t~ g{xt- a ) 



(20) 



where u = xt~ a is the similarity variable and a is the ex- 
ponent in the spreading law. Substituting this in Eq. I|19|) 
obtains both an exponent relation 



2a + (3 = 1 



(21) 



which must be satisfied for the similarity solution to hold, 
and an ordinary differential equation (ODE) for the sim- 
ilarity function 



fl!) 



+ 2{g') 2 + aug' + (3g = 0. 



(22) 



This is a non-linear second order ODE with boundary 
conditions g(0) = A (which is set by the drop volume) 
and g'{0) = (required by symmetry). 

A second exponent relation follows from the integrated 
conservation law. Volume conservation dictates that 17 oc 
J_ h 2 dx is constant. Inserting the similarity solution 

g 2 du is constant, 



in this shows that f2 
and therefore 



-26 



a = 2(3. 

Solving Eqs. (ED and 123 gives 

a = 2/5, /3 = l/5. 



(23) 



(24) 



Thus the spreading law L ~ t 2 ^ 5 of the preceeding sec- 
tion is recovered. The dependence on 17 and a/rj can be 
determined by dimensional analysis. 

To complete the discussion, the ODE for the similar- 
ity function g(u) can be solved. Inserting Eq. (|24|l into 
Eq. lO results in gg" + 2(g') 2 + 2ug'/h + g/5 = 0, with 
g(0) = A and <?'(0) = 0. Remarkably, this equation has 
an extremely simple closed form solution, 



g = A - u 2 /10. 



(25) 



Note that g — * for u — * no = V4JL4, so the spreading 
drop in the groove has a finite extent. Some examples of 
the shape for different values of A are shown in Fig. 0{a) • 
The basic prediction is that the height profile of the liquid 
surface for a droplet spreading in wedge is a parabola 
in the distance along the wedge. To be specific, from 
Eqs. (H3J), <Ei and (S5f. 



h(x,t) = h Q [1 - (x/x ) 2 }, (\x\<x ) 



(26) 



where the height ho ~ i -1 / 5 and the maximum extent of 
spreading xo ~ t 2 / b . 

For the case of spreading in a network of grooves, 
Eqs. I|19|) - (|22(l remain the same (apart from numeri- 
cal prefactors [lj]), but the integrated conservation law 
changes to 17 oc J °° h 2 2nrdr/l (where r should replace x 
in Eq. Ijl9|l ). Substituting the similarity solution shows 
that a = f3 must hold in this case. Combining this with 
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Eq. i(2"T|) gives a = (3 = 1/3, thus the previous scaling 
exponent is recovered. Unfortunately the corresponding 
ODE does not have a closed form solution. Some ex- 
amples of shapes obtained by numerical integration are 
shown in Fig. (2fb) . Again g — ► for u — + uq but in this 
case g has a singularity: g ~ (uo — u) 1 / 3 as u -> «o from 
below. 

Finally, the problem of a liquid spreading along a hy- 
drophilic strip shown in Fig. QJd) is discussed. The gov- 
erning partial differential equation has been obtained by 
Darhuber et al for this problem It is dh/dt ~ 

d[h 3 (dh/dx)]/dx where the prefactor can be found from 
Eq. (7) in Ref. 15] . Substituting the similarity trial so- 
lution Eq. J20I) in this yields an ODE g 3 g" + 3g 2 (g') 2 + 
aug' + f3g — with g(0) = A and g'(0) — being ini- 
tial conditions, and an exponent relation 2a + 3(3 = 1. 
It follows from f2 oc f_ wh dx that a = (3. Solving 
this together with the preceeding exponent relation gives 
a — (3 — 1/5, thus recovering the exponent of the pre- 
ceeding section. Numerical integration of the ODE gives 
results which resemble those for the case of spreading in 
a network of grooves, shown in Fig. |2Jb). Again g — > 
as u — > uo with g ~ (un — w) 1 / 3 as u — > uq from below 
(compare Fig. 3 of Ref. Ha) . 

In all these cases, the fact that g — > for u — > uq stands 
in contrast to the similarity solution for Tanner's law for 
which there is no point where g — ► Q. In that case 
one has to invoke an additional microscopic mechanism 
to account for the shape of the edge of the drop. 

The case of spreading from a reservoir can also be 
treated with variant of the above analysis; in fact the es- 
sential arguments are already given by Romero and Yost 
[l2|. and Darhuber et al [l^. In the case of a reservoir, h 
is constant at the reservoir edge which we define to be the 
point x = 0. In terms of the similarity solution Eq. (J2DJ, 
this forces (3 = 0. Combining this with the exponent re- 



lation in Eq. I|21|) or the analogous exponent relation for 
the strip problem shows a = 1/2 for all cases. Thus the 
Washburn-like spreading law is recovered for spreading 
from a reservoir, independent of whether spreading takes 
place in a wedge, in a network of grooves, or along a strip. 



IV. DISCUSSION 

The main results concern the kinetics of spreading in 
various geometries. New predictions are made for the 
scaling laws governing the rate at which a droplet spreads 
in a wedge or V-shaped groove, in a network of such 
grooves, and on a hydrophilic strip. These are estab- 
lished both by simple scaling arguments and by similarity 
solutions of the underlying partial differential equations. 
The asymptotic shapes of the spreading droplets have 
also been considered. 

Previous work on these problems has assumed the pres- 
ence of a reservoir which supplies liquid at constant pres- 
sure. This results in spreading laws which are essentially 
the same as the Washburn law for penetration of a liquid 
into a porous material, or into a capillary. The analysis 
here complements this previous work by considering the 
problems in the absence of a reservoir. This will apply 
in the late stages of wetting when the reservoir becomes 
exhausted, or if there is only a small amount of liquid 
present. 

These predictions could be tested by simulation or in 
experiments. The case of a droplet in a wedge would seem 
to be particularly simple, for instance a drop will spread 
in a right-angled corner provided the contact angle is less 
than 45°. The prediction that the scaling shape of the 
droplet should be a parabola should also be tested. 

I thank Alex Lips for discussions and encouragement. 
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